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CLAIMS 



What is claimed is: 

1 . A method of inverting a 4x4 source matrix, the method comprising: 
dividing the source matrix into four 2x2 sub-matrices A, B, C and D; 
calculating a plurality of sub-matrix products from the sub-matrices; 
calculating a determinant of the source matrix dS to form a matrix determinant 

residue rd of the source matrix as rd=l/dS; 

forming a partial, inverse sub-matrix of each sub-matrix using one or more of the 
matrix products and a determinant of each sub-matrix; and 

calculating an inverse of each sub-matrix iA, iB, iC, and iD, utilizing each partial, 
inverse sub-matrix and determinant residue rd, such that an inverse of the source matrix iS 
is formed. 

2. The method of claim 1, wherein dividing the source matrix S into the four 
2x2 sub-matrices A, B, C and D is performed according to the following rule: 



to enable storage of each sub-matrix within a pair of SIMD registers. 

3 . The method of claim 1 , wherein calculating the plurality of sub-matrix 
products further comprises: 

calculating an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 



wherein the adj function refers to an adjoint matrix operation and the dot symbol • 
refers to a matrix multiplication operation; and 




DC = adj(D) • C 
A5 = adj(A)«B 
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calculating a final sub-matrix product for each of the intermediate sub-matrix 
products by computing the following equations: 
BDC = B*DC 
DBA = D«adj(AB) 
ACD = A»adj(DC) 
CAB = C •AB. 

4 . The method of claim 1 , wherein calculating the matrix determinant residue 
further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD; 
calculating a trace value by computing a following equation: 

t = trace {AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; and 

calculating a determinant of the source matrix dS by computing a following 
equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation. 

5 . The method of claim 1 , wherein forming partial -inverse sub-matrices further 
comprises: 

performing matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, 
B*dC and A*dD; and 

computing a partial inverse for each sub-matrix according to the following matrix 
scaling equations: 

pA = A*dD-BDC 
pB = C*dB- DBA 
pC=B*dC-ACD 
pD = D*dA- CAB, 
wherein pA, pB, pC, and pD reference partial, inverse sub-matrices, and the 
symbol * refers to a matrix scaling by a scalar operation. 
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6 . The method of claim 1, wherein calculating an inverse of each sub-matrix 
further comprises: 

calculating an adjoint value of each partial, inverse sub-matrix pA, pB, pC, and pD, 
according to the following rules: 
iA = adj(pA), 
iB = ad}(pB), 
/C=adj(pC), 
iD = ad)( P D), 

wherein the adj() function refers to the adjoint matrix operation; 

calculating a final sub-matrix inverse value according to the following equations: 

iA = iA * rd 

iB = iB * rd 

iC = iC * rd 

iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
forming the inverse source matrix iS according to the following rule: 



7 . A method compri sing: 

dividing a source matrix into four 2x2 sub-matrices, A, B, C and D; 
calculating one or more intermediate sub-matrix products from one or more of the 
sub-matrices; 

calculating a determinant of the source matrix to form a determinant residue rd 
utilizing the intermediate sub-matrix products; 

scaling a determinant of each sub-matrix and the intermediate sub-matrix products 
using determinant residue rd to form final sub-matrix products; 

forming a partial inverse sub-matrix pA, pB, pC and pD for each sub-matrix using 
the scaled sub-matrix determinants and the final sub-matrix products; and 

calculating an inverse of each sub-matrix iA, iB, iC and iD, utilizing each partial 
inverse sub-matrix to form an inverse source matrix iS. 
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8 . The method of claim 7, wherein calculating the matrix determinant residue 
further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD; 
calculating a trace value by computing a following equation: 

t = trace (AB»DC); 
wherein a dot symbol • refers to a matrix multiplication operation; 

calculating a determinant of the source matrix dS by computing a following 
equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
calculating the determinant residue rd according to the following rule: 

rd = 1/dS. 

9 . The method of claim 7, wherein scaling by the determinant residue further 
comprises: 

multiplying each determinant by the determinant residue rd according to the 
following rules: 

dA = dA * rd 
dB = dB * rd 
dC = dC * rd 
dD = dD * rd; 

multiplying each intermediate sub-matrix product AB and DC by the determinant 
residue rd, according to the following equations: 
DC=DC*rd 
AB = AB * rd; and 

calculating a final sub-matrix product for each of the intermediate matrix products 
by computing the following equations: 
BDC = B*DC 
DBA = D«adj(AB) 
ACD = A»adj(f»c) 
CAB = C •AB. 



28 



Ref. No! 42390P13115 



Express Mail No.: EL651893785US 



10. The method of claim 7, wherein calculating an inverse of each sub-matrix 
further comprises: 

generating an adjoint of each partial, inverse sub-matrix by computing the following 



iA = adj(pA) 
iB = adj(pfi) 
iC = adj(/?C) 
iD = adj(pD); and 
forming the inverse source matrix iS according to the following rule: 



11. A computer readable storage medium including program instructions that 
direct a computer to function in a specified manner when executed by a processor, the 
program instructions comprising: 

dividing the source matrix into four 2x2 sub-matrices A, B, C and D; 

calculating a plurality of sub-matrix products from the sub-matrices; 

calculating a determinant of the source matrix dS to form a matrix determinant 
residue rd of the source matrix as rd=l/dS; 

forming a partial, inverse sub-matrix of each sub-matrix using one or more of the 
matrix products and a determinant of each sub-matrix; and 

calculating an inverse of each sub-matrix iA, iB, iC, and iD, utilizing each partial, 
inverse sub-matrix and determinant residue rd, such that an inverse of the source matrix iS 
is formed. 

1 2 . The computer readable storage medium of claim 1 1 , wherein dividing the 
source matrix S into the four 2x2 sub-matrices A, B, C and D is performed according to the 
following rule: 



equations: 





to enable storage of each sub-matrix within a pair of SIMD registers. 
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1 3 . The computer readable storage medium of claim 1 1 , wherein calculating the 
plurality of sub-matrix products further comprises: 

calculating an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 



AB = adj( A) •£ 

wherein the adj() function refers to an adjoint matrix operation and the dot symbol • 
refers to a matrix multiplication operation; and 

calculating a final sub-matrix product for each of the intermediate sub-matrix 
products by computing the following equations: 
BDC = B»DC 



CAB = C •AB. 

1 4 . The computer readable storage medium of claim 1 1 , wherein calculating the 
matrix determinant residue further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD; 
calculating a trace value by computing a following equation: 

t = trace {AB • DC); 
wherein a dot symbol • refers to a matrix multiplication operation; and 
calculating a determinant of the source matrix dS by computing a following 





equation: 



dS = dA*dD + dB*dC - t 



wherein the symbol * refers to a scalar multiplication operation. 
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1 5 . The computer readable storage medium of claim 1 1 , wherein forming 
partial-inverse sub-matrices further comprises: 

performing matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, 
B*dC and A*dD; and 

computing a partial inverse for each sub-matrix according to the following matrix 
scaling equations: 



wherein pA, pB, pC, and pD reference partial, inverse sub-matrices, and the 
symbol * refers to a matrix scaling by a scalar operation. 

1 6 . The computer readable storage medium of claim 1 1 , wherein calculating an 
inverse of each sub-matrix further comprises: 

calculating an adjoint value of each partial, inverse sub-matrix pA, pB, pC, and pD, 
according to the following rules: 
iA = ad}(pA), 
iB = adj(pB), 
iC = adj(pC), 
iD = adj(pD), 

wherein the adj() function refers to the adjoint matrix operation; 

calculating a final sub-matrix inverse value according to the following equations: 

iA = iA * rd 

iB = iB * rd 

iC = iC * rd 

iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
forming the inverse source matrix iS according to the following rule: 



pA. 

P B: 

pC 
pD 



A*dD-BDC 



C*dB- DBA 



B*dC-ACD 



D*dA- CAB, 
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1 7 . The computer readable storage medium including program instructions that 
direct a computer to function in a specified manner when executed by a processor, the 
program instructions comprising: 

dividing a source matrix into four 2x2 sub-matrices, A, B, C and D; 
calculating one or more intermediate sub-matrix products from one or more of the 
sub-matrices; 

calculating a determinant of the source matrix dS to form a determinant residue rd 
of the source matrix utilizing the intermediate sub-matrix products and the sub-matrix 
determinants; 

scaling a determinant of each sub-matrix and the intermediate sub-matrix products 
using determinant residue rd to form final sub-matrix products; 

forming a partial inverse sub-matrix pA, pB, pC and pD for each sub-matrix using 
the scaled sub-matrix determinants and the final sub-matrix products; and 

calculating an inverse of each sub-matrix iA, iB, iC and iD, utilizing each partial 
inverse sub-matrix to form an inverse source matrix iS. 

1 8 . The computer readable storage medium of claim 17, wherein calculating the 
matrix determinant residue further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD; 
calculating a trace value by computing a following equation: 

t = trace (AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; 

calculating a determinant of the source matrix dS by computing a following 
equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
calculating the determinant residue rd according to the following rule: 

rd = 1/dS. 
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19. The computer readable storage medium of claim 17, wherein scaling by the 
determinant residue further comprises: 

multiplying each determinant by the determinant residue rd according to the 
following rules: 

dA = dA * rd 
dB = dB * rd 
dC = dC * rd 
dD = dD * rd; 

multiplying each intermediate sub-matrix product by the determinant residue rd, 
according to the following equations: 
DC = DC*rd 
AB = AB*rd; and 

calculating a final sub-matrix product for each of the intermediate matrix products 
by computing the following equations: 
BDC = B*DC 
DBA = D«adj(l5) 
ACD = A«adj(Dc) 
CAB = C *AB. 

20. The computer readable storage medium of claim 17, wherein calculating an 
inverse of each sub-matrix further comprises: 

generating an adjoint of each partial, inverse sub-matrix by computing the following 
equations: 

iA = adj(pA) 
iB = adj(pB) 
iC = adj(pC) 
iD = adj(pD); and 
forming the inverse source matrix iS according to the following rule: 

{iC iDJ. 
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21. An apparatus, comprising: 

a processor having circuitry to execute instructions; 

a plurality of SIMD data storage devices coupled to the processor, the SIMD data 
storage registers to pairs of floating point vectors during matrix calculation; 

a storage device coupled to the processor, having sequences of instructions stored 
therein, which when executed by the processor cause the processor to: 

divide the source matrix into four 2x2 sub-matrices A, B, C and D; 

calculate a plurality of sub-matrix products from the sub-matrices; 

calculate a determinant of the source matrix dS to form a determinant residue 
rd of the source matrix as rd=l/dS; 

form a partial, inverse sub-matrix of each sub-matrix using one or more of 
the matrix products and the determinant of each sub-matrix; and 

calculate an inverse of each sub-matrix iA, iB, iC, and iD, utilizing each 
partial, inverse sub-matrix and determinant residue rd, such that an inverse of the source 
matrix iS is formed. 

22. The apparatus of claim 21 , wherein the instruction to calculate the plurality 
of sub-matrix products further causes the processor to: 

calculate an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 

DC = adj(z>) • C 

AB = adj(A) • B 

wherein the adj() function refers to an adjoint matrix operation and the dot symbol • 
refers to a matrix multiplication operation; and 

calculate a final sub-matrix product for each of the intermediate sub-matrix products 
by computing the following equations: 
BDC = B*DC 
DBA = D*ad^ABj 
ACD = A«adj(Dc) 
CAB = C •AB. 
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23. The apparatus of claim 21, wherein the instruction to calculate the matrix 
determinant residue further causes the processor to: 

compute a determinant of each sub-matrix dA, dB, dC and dD; 
calculate a trace value by computing a following equation: 

t = trace (AB •DC); 
wherein a dot symbol • refers to a matrix multiplication operation; and 
calculate a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation. 

24. The apparatus of claim 21, wherein the instruction to perform matrix scaling 
further causes the processor to: 

perform matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, B*dC 
and A *DdD; 

compute a partial inverse for each sub-matrix according to the following matrix 
scaling equations: 

P A = A*dD-BDC 
pB — C*dB- DBA 
pC= B*dC-ACD 
pD = D* dA - CAB, 
wherein pA, pB, pC, and pD reference partial, inverse sub-matrices and the 
symbol * refers to a matrix scaling by a scalar operation. 

25 . The apparatus of claim 21 , wherein the instruction to calculate an inverse of 
each sub-matrix further causes the processor to: 

calculate an adjoint value of each partial, inverse sub-matrix pA, pB, pC, andpD, 
according to the following rules: 
iA = adj(pA), 
iB = adj(pB), 
*'C = adj(pC), 
iD = adj(pD), 

wherein the adj() function refers to the adjoint matrix operation; 
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calculate a final sub-matrix inverse value according to the following equations: 
iA = iA * rd 
iB = iB * rd 
iC = iC * rd 
iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
form the inverse source matrix iS according to the following rule: 

is J" iB ) 

[iC iD J . 

26. An apparatus, comprising: 

a processor having circuitry to execute instructions; 

a plurality of SIMD data storage devices coupled to the processor, the SIMD data 
storage registers to pairs of floating point vectors during matrix calculation; 

a storage device coupled to the processor, having sequences of instructions stored 
therein, which when executed by the processor cause the processor to: 

divide a source matrix into four 2x2 sub-matrices, A, B, C and D; 
calculate one or more intermediate sub-matrix products from each of the 

sub-matrices, 

calculate a source matrix dS to form a determinant residue rd utilizing the 
intermediate sub-matrix products, 

scale a determinant of each sub-matrix and the intermediate sub-matrix 
products using determinant residue rd to form final sub-matrix products, 

form a partial inverse sub-matrix pA, pB, pC and pD for each sub-matrix 
using the scaled sub-matrix determinants and the final sub-matrix products, and 

calculate an inverse of each sub-matrix iA, iB, iC and iD, utilizing each 
partial inverse sub-matrix to form an inverse source matrix iS. 
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27 . The apparatus of claim 26, wherein the instruction to calculate the source 
matrix determinant residue further causes the processor to: 

compute a determinant of each sub-matrix dA, dB, dC and dD; 
calculate a trace value by computing a following equation: 

t = trace (AB»DC) 
wherein a dot symbol • refers to a matrix multiplication operation; 
calculate a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
calculate the determinant residue rd according to the following rule: 

rd = 1/dS. 

28. The apparatus of claim 26, wherein the instruction to scale by the 
determinant residue further causes the processor to: 

multiply each determinant by the determinant residue rd according to the following 

rules: 

dA = dA * rd 
dB = dB * rd 
dC = dC * rd 
dD = dD * rd; 

multiply each intermediate sub-matrix product AB and DC by the determinant 
residue rd, according to the following equations: 
DC = DC*rd 
AB = AB * rd; and 

calculate a final sub-matrix product for each of the intermediate matrix products by 
computing the following equations: 

BDC = B*DC 

DBA = D«adj(AB) 

ACD = A»adj(z)c) 

CAB = C • AB. 
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29. The apparatus of claim 26, wherein the instruction to calculate an inverse of 
each sub-matrix further causes the processor to: 

generate an adjoint of each partial, inverse sub-matrix by computing the following 



iA = adj(jpA) 
iB = adj(pB) 
iC = adj(pC) 
iD = adj(pD); and 
form the inverse source matrix iS according to the following rule: 



equations: 
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